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The quantum ferromagnetic transition at zero temperature in disordered itinerant electron systems 
is considered. Nonmagnetic quenched disorder leads to diffusive electron dynamics that induces 
an effective long-range interaction between the spin or order parameter fluctuations of the form 
r 2_2d , with d the spatial dimension. This leads to unusual scaling behavior at the quantum critical 
point, which is determined exactly. In three-dimensional systems the quantum critical exponents 
are substantially different from their finite temperature counterparts, a difference that should be 
easily observable. Experiments to check these predictions are proposed. 



I. INTRODUCTION 

The general problem of describing quantum phase 
transitions is a subject of great current interest. These 
transitions occur at zero temperature as a function of 
some nonthermal control parameter, and the relevant 
fluctuations are of quantum nature rather than of ther- 
mal origin. Early work in this field established that if 
the quantum phase transition has a classical analogue 
at finite temperature, then in the physical dimension, 
d = 3, the former tends to have a simpler critical behav- 
ior than the latter. In particular, one of the most obvi- 
ous examples of a quantum phase transition, namely the 
ferromagnet-to-paramagnet transition of itinerant elec- 
trons at zero temperature, T — 0, as a function of the ex- 
change interaction between the electron spins,.-svas found 
to have mean-field critical behavior in d = 3u The rea- 
son for this tendency is that the coupling between stat- 
ics and dynamics that is inherent to quantum statistics 
problems effectively increases the dimensionality of the 
system from d to d + z, with z the dynamical critical 
exponent. In the case of clean itinerant quantum ferro- 
magnets, z — 3 in mean-field theory. This appears to 
reduce the upper critical dimension d+, which is the di- 
mension above which mean-field theory yields the exact 
critical behavior, from <i+ =„4 in the classical case to 
d£ = 1 in the quantum case,El However, this conclusion 
has recently been disputedpfl and the critical behavior of 
clean itinerant ferromagnets in low-dimensional systems 
is currently under renewed investigation^ 

It has been known for some time that if one adds 
quenched, nonmagnetic disorder to the system, then the 
critical behavior at the quantum phase transition must be. 
substantially modifiedfl contrary to earlier suggestions .u 
It is known that the correlation length exponent, v, must 
satisfy the inequality v > 2/d in systems with quenched 



disorder The physical origin of this constraint is the 
requirement that disorder induced fluctuations in the lo- 
cation of the critical point in parameter space must be 
small compared to the distance from the critical point in 
order for the phase transition to be sharp. Any mean- 
field, and any standard Gaussian theory yields v = 1/2, 
which is incompatible with this lower bound in any di- 
mension d < 4. Technically, this implies that the disorder 
is a relevant perturbation with respect to the mean-field 
or Gaussian fixed point. 

From a perspective that is entirely focused on the sta- 
tistical mechanics aspects of the phase transition problem 
(and that therefore does not take into account from the 
beginning the aspects that have to do with the under- 
lying disordered electron system), it would be tempting 
to model this disordered quantum phase transition by 
making the mass term in the effective field theory for the 
clean case (i.e. the coefficient of the term quadratic in 
the order parameter field in the Landau-Ginzburg- Wilson 
(LGW) functional) a random variable in order to describe 
the fluctuations in the location of the critical point. For 
the quantum paramagnet-to-ferromagnet transition un- 
der consideration one can readily convince oneself that 
for d < 4 the resulting disorder term is relevant in the 
renormalization group sense with respect to the clean 
Gaussian fixed point. Presumably, the presence of this 
relevant operator leads to a new critical fixed point with 
a correlation length exponent v > 2/d. 

It turns out, however, that such a model is not a tech- 
nically appropriate description of the quantum ferromag- 
netic transition in a system of disordered itinerant elec- 
trons. In order to explain this important point, let us 
anticipate a number of results that will be discussed in 
detail in Sec. || below. If one attempts to derive the 
LGW functional mentioned above, then the coefficient of 
the random mass term is a correlation function whose 
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expansion in the random potential is very singular. The 
reason is that the dynamics of a quantum particle in a dis- 
ordered environment in the limit of small wavenumbers, 
k, and low frequencies, u), are qualitatively different from 
those in a clean system. In the former case the motion 
is diffusive, while in the latter it is ballistic. Technically, 
the limit fc, uj — > does not commute with the clean limit. 
To overcome this problem it is advantageous to not ex- 
pand in powers of the random potential, but to expand 
instead in the fluctuations of the coefficients of the LGW 
functional at fixed disorder. 

This procedure, by construction, automatically resums 
the most obvious singularities in the disordered itinerant 
electron problem. Although it will turn out that the coef- 
ficients in the resulting LGW functional are still singular, 
these remaining singularities can be handled mathemati- 
cally as we will see below. These singularities arise. due to 
what are often called weak localization effectsE'El Their 
physical origin lies in the fact that the dynamics of dif- 
fusive electrons are intrinsically long-ranged. Via mode- 
coupling effects, this feature couples back to the ferro- 
magnetic transition problem. As we show in this paper, 
the net effect is that the interactions between the spin 
fluctuations in the LGW theory are of long range. The 
critical behavior described by the resulting nonlocal field 
theory can be determined exactly for all d > 2, and satis- 
fies v > 2/d as required. The critical exponents obtained 
from this theory are ci-dependent for all d < 6. In d = 3, 
they are substantially different from either the mean-field 
exponents, or from those for a classical Heisenberg model, 
which has striking observable consequences. 

The outline of this paper is as follows. In Sec. |l| we first 
define an itinerant disordered electron model, and then 
discuss how to derive an order parameter description for a 
ferromagnetic phase transition starting from a fermionic 
field theory. We also discuss in detail the coefficients in 
the LGW functional and show that they have long-range 
properties. In Sec. [II the critical behavior is determined 
exactly for dimensions d > 2. For d > 6, mean- field ex- 
ponents are obtained, while for 2 < d < 6, d-dependent 
exponents are found. In Sec. IV the results of this paper 
are reviewed. First we summarize the theoretical aspects 
of our results, and then we point out some of the ex- 
perimental consequences we expect. In Appendix A we 
calculate the wavenumber dependent spin susceptibibil- 
ity for disordered interacting fermions, and in Appendix 
B we discuss the logarithmic corrections to scaling that 
exist at the critical dimensions d^ = 6 and d^ = 4, as 
well as for all 2 < d < 4. 



with the spin density fluctuation as the order parameter. 
In the last part of this section we derive and discuss the 
coefficients in this LGW functional. As mentioned in the 
introduction, the crucial point is that the interactions 
in the effective LGW theory are long-ranged due to the 
diffusive dynamics of the electrons in a disordered metal. 



A. The Model 

The partition function of any fermionic system can be 
written in the form,Ll 



Z = J Dip Dip exp(S[i>,ip]) 



(2.1a) 



where the functional integration measure is defined with 
respect to Grassmannian (i.e., anticommuting) fields tp 
and tp, and S is the action, 

S= f dr J dx^(x,r)^-^(x,r)- f dr H{t) . 

(2.1b) 

Here x denotes positions and r imaginary time, H(t) 
is the Hamiltonian in imaginary time representation, 
j3 = 1/T is the inverse temperature, i = 1,2 denotes 
spin labels, and summation over repeated indices is im- 
plied. Throughout this paper we use units such that 
ks = h = e 2 = 1. Our starting model is a fluid of in- 
teracting electrons moving in a static random potential 
«(x), 



H(t) 



rfx 



1 

2m 



V^(x,t) • V#(x, r) 



+[w(x) - n] ^(x, t) V> l (x, T ) 



+ - j dxdy u(x-y) tp 1 (x, t) tp 1 (y, t) tp 1 (y, r) 

x^(x,r) . (2.2a) 

Here m is the electron mass, [i is the chemical potential, 
and u(x — y) is the electron-electron interaction poten- 
tial. We assume that the random potential u(x) is delta- 
correlated and obeys a Gaussian distribution P[v(x)] 
with second moment 



{v{x)v(y)} di 



1 



2TTN F T e 



5(x-y) , 



(2.2b) 



II. MODEL, AND THEORETICAL FRAMEWORK 

In the first part of this section we define a simple 
model for interacting electrons in a disordered environ- 
ment. Starting with this general fermionic field theory 
we then derive an LGW or order parameter description, 



where 



V ■ -}dis 



D[v] P[v] (. . .) 



(2.2c) 



denotes the disorder average, Np is the bare density of 
states per spin at the Fermi level, and r e ; is the bare 
electron elastic mean-free time. More realistic models 
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to describe itinerant electron magnetism including, e.g., 
band structure, can be considered along the same lines. 
The salient points of our results, however, are due to 
long-wavelength effects and hence do not depend on mi- 
croscopic details like the band structure. For our pur- 
poses it therefore is sufficient to study the model defined 
in Eqs. 

For describing magnetism, it is convenient and stan- 
dard practice to break the interaction part of the action 
5, which we denote by Si n t, into spin-singlet and spin 

(s t) 

triplet contributions, S^ n l . For simplicity, we assume 
that the interactions are short-ranged in both of these 
channels. In a metallic system this is justified due to 
screening, and an effective model with a short-ranged in- 



teraction in both S^l and S£i can be derived starting 
from a bare Coulomb interaction.EHl In order for this as- 
sumption to remain valid, our discussion applies only to 
cases in which the disorder is weak enough for the sys- 
tem to remain far from any metal-insulator transition 
that might be present in the phase, diagram and would 
lead to a breakdown of screening.!!! The spin-triplet in- 
teraction S^t describes interactions between spin density 
fluctuations. This is what causes ferromagnetism, and it 
therefore makes sense to consider this part of the action 
separately. We thus write 



(t) 



S = So + S 



with 



5' 



(t) 
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J dxdr n s (x, r) • n s (x,' 



(2.3) 



(2.4a) 



where n s is the electron spin density vector with compo- 
nents, 



n:(x,r)^f(x,T)<y(x,T) 



(2.4b) 



Here the cr° are the Pauli matrices, and T t is the spin 
triplet interaction amplitu de th at is related to the inter- 
action potential u in Eq. (2.2a) via 



r, = 



<ix u(x) 



(2.4c) 



For simplicity we have assumed a point-like interaction 
so that r f is simply a number. A generalization to a 
more realistic short- ran ge interaction would be straight- 
forward. Sq in Eq. (2.3) contains all other contributions 
to the action. It reads explicitly, 



S () 



dr / d: 



d 



2m 



V^(x,t) ■ W/(x,t) 



with n c the electron charge or number density, 

n c (x,r) =#(x,T)V i (x,-r) , (2.5b) 
and T s the spin-singlet interaction amplitude. 

B. Order parameter field theory 

Continuous thermal phase transitions are usually de- 
scribed by deriving an LGW theory, i.e. an effective 
field theory for the long wa\jelength order parameter fluc- 
tuations or critical modes. t2l The physical idea behind 
this approach is that these fluctuations, which are slowly 
varying in space, determine the behavior near the critical 
point. The same philosophy has been applied to quantum 
phase transitions, with the only difference being that the 
critical modes are now slowly varying in both space and 
time. We will use this approach here, motivated in pari, 
by previous work on clean itinerant electronic systems .EJ 
We mention, however, that in general one should worry 
about both the critical modes, and all other slow or soft 
modes, even if these other soft modes are not 'critical' in 
the sense that they change their character at the phase 
transition. While this concern is not confined to quan- 
tum phase transitions, we will argue below that for these 
it poses a more serious problem than for thermal phase 
transitions since at T = there are more soft modes than 
at finite temperature. We will see that in the present 
problem such additional modes are indeed present and 
lead to complications within the framework of an LGW 
theory. For the problem under consideration, however, 
these complications can be overcome. 

The techniques for deriving an order para mete r field 
theory, starting with Eqs. (|2.3D through (|2.5|), are 



standard.! We decouple the four fcrmion term in S£t by 
introducing a classical vector field M(x, r) whose average 
is proportional to the magnetization m, and performing 
a Hubbard-Stratonovich transformation. All degrees of 
freedom other than M are then integrated out. This 
procedure in particular integrates out the soft diffusive 
modes or 'diffusotis' that are inherent to a disordered 
fermion system.LM3 These are the additional soft modes 
mentioned above. We obtain the partition function Z in 
the form, 



(*) 



Z = e - Fo/T / D[M] exp(-$[M]) 



(2.6a) 



with Fq the noncritical part of the free energy. The LGW 
functional <i> reads, 



$[M] = y J dx M(x) ■ M(a 



[u(x)-/x]V i (x,r) ^(x,t) 



In yexp 



If 
2 



-T t J dx M(x) ■ n s (x) 



fP r where we nave adopted a tour-vector 

^ drj dxn c (x,r)n c (x,r) , (2.5a) (X)T)j and J & = J dx j/V. Here 



(2.6b) 



So 



where we have adopted a four-vector notation with x — 

denotes an 



)s 
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average taken with the action Sq. A formal expansion of 
$ in powers of M takes the form, 



$[M] = i J dxidx 2 M a {xi) 



— 6(xi -x 2 ) 
1 t 



M b {x 2 ) 



+ 



1 

4! 



^ J dxidx 2 dx 3 Xall( x l, x 2,X 3 ) 

xM a ( Xl )M b (x 2 )M c (x 3 ) 

dxi dx 2 dx 3 dxi Xabld( x ~>--> x ^ x 3, 

x M a ( Xl ) M b {x 2 ) M c {x 3 ) M d { Xi ) 
+ 0(M 5 ) , 



(2.7a) 

where we ha ve scaled M with T^ 1 . The coefhcients x^ 



in Eq. ( 2.7a ) are correlation functions for a system with 
a particular realization of the disorder (i.e., they are not 
translationally invariant). They are defined as, 



(2.7b) 



where the superscript 'c' denotes a cumulant or con- 
nected correlation function. For our simple model, the 
reference ensemble with action So, whose correlation 
functions are the x j consists of free electrons with dis- 
order and a short-ra nged spin-singlet model interaction. 
As mentioned in Sec. II A above, the model can be made 
more realistic, if desired, by e.g. including a realistic 
band structure. In order to do so, one would simply re- 
place the x above with correlation functions for band 
electrons. 



Equations (2.6) and (2.7) define an order parameter 
field theory for the paramagnet-to-ferromagnet phase 
transition in a disordered itinerant electronic system. 
The disorder has not been averaged over yet, and so the 
rt-point correlation functions in Eqs. (2.7) depend explic- 
itly on the particular realization of the randomness in the 
system. To proceed, we now formally carry out the disor- 
der average of the free energy, using replicas and a cumu- 
lant expansion, and keeping—terms up to 0(M 4 ) in the 
replicated LGW functional^ With a, /3, . . . = 1, . . . , n 
(n — ► 0) the replica labels, we obtain the following LGW 
functional for the a th copy of the replicated system, 



$ Q [M] 



5>f[M] 

i= 2 



= - j d Xl dx 2 X%(x u x 2 )M:( Xi )M«(x 2 ) 



+ 
1 

4! 



dxi dx 2 dx 3 X x ab ' c (x!,x 2 , x 3 ) 

xM^x^M^x^M?^) 
dxidx 2 dx 3 dx A X^ P (x 1 ,x 2 ,x 3 ,x 4: ) 

xM:( Xi ) M?(x 2 ) MfOca) M$(X4) 

+0(M 5 ) , (2.8a) 



where <f>" denotes the contribution of order M , and the 
coefficients are given in terms of disorder averaged corre- 
lation functions, 



Xf h \x x ,x % ) 



~^-5{xi -X 2 ) - {X [ ab\xi,x 2 )} dl 
1 t 



(2.8b) 



X^ b \(x u x 2 ,x 3 ) = {xj^l (x!,x 2 , x a )}di s , (2.8c) 

X^l^(xi,X 2 ,X 3 ,X4,) = S a p{Xai cd (xi,X 2 ,X 3 ,X i )} dis 

+ 3{x^(x 1 ,x 2 )x^(x 3 ,x i )}' dls . (2.8d) 

Again, the superscript 'c' means that only the connected 
part, this time with respect to the disorder average, of 
the correlation functions {. . .} dis should be considered. 
Notice that we have separated the quantum mechani- 
cal and disorder averages, and therefore must deal with 
two different kinds of cumulants: In Eqs. ( p77|) , we have 
{ab) c = (ab) — (a) (b), etc. while in Eqs. Q2.8| ) we have 
{ ab }dis = {ab}di S - {a}dis{b}dis, etc. 



C. The coefficients 



As it stands, Eq. (2.8a) is just a formal nonlocal field 
theory that is not very useful. Normally, one would 
proceed by localizing the individual terms in Eq. (2.8a) 



about a single point in space and time, and expanding the 
correlation functions in powers of gradients. Due to what 
are often called weak localization effects jlil this is not an 
option for the current case of a disordered metal, even 
if we stay far away from any metal-insulator transition. 
Structurally it is easy to show, using either many-body 
perturbation theory or field theoretic methods, that the 
disorder averaged n-point correlation function, {x }disj 
in wavenumber space has a weak localization correction 
of a form that can be written schematically as, 



f 



k>q uj n >0„ 



(k 2 



(2.9) 



with V the system volume, k — |k|, u) n and f2„ Mat- 
subara frequencies, and (q, fi n ) an infrared cutoff. Al- 
though, strictly speaking, <5x depends on I — I external 
momenta and frequencies, we have schematically repre- 
sented these by a single 'typical' momentum-frequency 
(q, fl n ). Since we will only be interested in the scale 
dimension of {x }dis, this is sufficient. In terms of di- 
agrams for the field theory reviewed in Ref. ||, the dom- 
inant contribution to {x^}dis{l, = 0) f° r 1 ~~ * is 
shown in Fig. [II For what follows it is very important to 
notice that 5x^'(Q^n — 0) is a nonanalytic function of 
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/ (n-6) 
times 



FIG. 1. Diagrammatic structure of the leading IR singular 
contribution to the correlation function x (q> = 0). The 
straight lines denote the diffusive propagators of the field the- 
ory reviewed in Ref. ^, and the dotted lines denote external 
legs. 



To illustrate this point in detail, let us consider the 
two-point correlation function, xi&\ nl Eq. ( 2.8b| ). Be- 
cause of translational invariance on average in space and 
time, as well as rotational invariance in space, we can 
write, 



( 2) 

{Xab( X l' X mdis = Sa,bXs(\Xl-X2\,n-T2) , (2.10a) 

and we define the Fourier transform 

Xs (cl, ti n ) = J d(xi - x 2 ) d(n - r 2 ) e"^'^-^) 

XetfWn-Ta) Xs ( Xl - X2 ) . (2.10b) 

Xs is the disorder averaged spin susceptibility of the ref- 
erence system whose action is given by So- Particle 
number conservation implies that at small frequency and 
wavenumber it has a diffusive structure, 



X s (q, ti n ) = xo(q) 



Dq 2 



\tir, 



Dq 2 



(2.11a) 



where D is the spin diffusion coefficient of the refer- 
ence ensemble, and Xo(l) is the static spin susceptibility. 
In a system with a conserved order parameter, the fre- 
quency must be taken to zero before the wavenumber in 
order to reach criticality, so in the critical region we have 
<< Dq 2 . Note that it is our restriction to nonmag- 
netic disorder that ensures a conserved order parameter. 
We will come back to this point below. In the critical 
limit, we can thus expand, 



Xs(cL,ti n ) = X0 (q) [1- \ti n \/Dq 2 + 



(2.11b) 



The static spin susceptibility, xo (q) > has been discussed 
in some detail before.Ej For any system with a nonvan- 
ishing interaction amplitude in the spin-triplet channel, 
r t =/= 0, there is a diffusive or weak localization correction 
to the bare susceptibility. Although our bare reference 



ensemble does not contain such an interaction amplitude, 
perturbation theory will generate one as long as there is a 
nonvanishing interaction amplitude r s in the spin-singlet 
channel.EJ Effectively, we therefore need to use the spin 
suscep tibility for a system with a nonvanishing T t in Eq. 
([2.111: ) . The explicit calculation is given in Appendix A, 



and the result is, 



Xo(q -> 0) = c - c d ^ 2 q 



d-2 



c 2 q 



(2.12a) 



where the Cj are positive constants. The nontrivial, and 
for our purposes most interesting, contribution in Eq. 
( |2.12a ) is the nonanalytic term ~ q d ~ 2 . Its existence 
implies tha t the standard gradient expansion mentioned 
above Eq. (2J3) does not exist. The physical interpre- 
tation of this term is that, effectively, there is a long- 
range interaction between the order parameter fluctua- 
tions, which in real space takes the form r ~ 2d+2 . This 
is a phenomenon that is special to zero temperature. At 
finite temperature, when one has to perform a frequency 
sum rather than a frequency integral to calculate the cor- 
relation funcion, the nonanalytic term is replaced by a 
term of the schematic structure, 



q 



(q 2 + T) (d_2)/2 



(2.12b) 



so for T > an analytic expansion about q = exists, 
and the standard local LGW functional is obtained. 

Rota- 



Next we consider the cubic term in Eq. (2 



tional symmetry in spin space allows us to write $3 as, 

$3 [M] = ^ J dxidx 2 dx 3 u 3 (x 1 ,x 2 ,x 3 ) 

xM a (n) • (M Q (a;2) x M"^)) , (2.13a) 



with 



u 3 (x 1: x 2 ,x 3 ) = - e abc {x y ab ' c {xi,x 2 ,x 3 )} dls , (2.13b) 
wher e e a bc is the Levi-Civita tensor. According to Eq. 



u 3 (q^0,ti n = 0)=u ( t 4) 



,(0) 



(2.13c) 



with u 3 d ^ and uf^ finite numbers. Note that this term 
is divergent for q — > for d < 4. 

The quartic term in Eq. ( [2.8a ) is written as a sum of 
two terms, 



$?[M] =< (1) [M]+$^ ZJ [M] 



a (2), 



(2.14a) 



with 



$4 W [M] 



dxi . . . dx A {xibL( x i' ■ ■ ■ ' x i)}di 



24 

xM^xr) M?(x 2 ) M?(x 3 ) M%{ Xi ) 



(2.14b) 
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and 



$4 ^ [M] = i J dxi dx 2 dx 3 dx 4 

X {Xab (XI , X 2 ) Xcf ( X 3 > a; 4)}dj 

xilCOn) M£(x 2 ) M?(x 3 ) M?( Xi ) 



(2.14c) 



For our purposes we only need to know the degrees of di- 
vergence of the coefficients in Eqs. (2.14). We can there- 
fore formally consider the local, static limit of Eq. ( 2.141 ) 
(even though it does not necessarily exist), and write 



<(i) 



[M] 



T 
2W 



"4 



E E 

qi,q2,q3 ni,n2,n3 



x(M Q (q 1 , W „ 1 )-M Q (q 2 , W „ 2 )) 
x (M Q (q3,^ n3 ) • M Q (- qi - q 2 - q 3 , -w ni 

-u n2 - oj„ 3 )) , (2.15a) 



where 



and 

W4 = 



r lim {xiaL}di S (qi, Q2, qs; w ni 



According to Eq. fcjh 



u i{q — * 0, Q„ = 0) = u 



- rv> - „( d " 6 ) „rf-6 



,(o) 



(2.15b) 



n 3 =0). 

(2.15c) 



(2.15d) 



so W4 diverges in this limit for d < 6, and is finite for 
d > 6. W ith a more accurate representation for m than 
Eq. (2.9) one also finds a term ~ q d ~ 4 , but this will be 
of no relevance for what follows. 

ft" ( 2 



We now consider <I> 4 given by Eq. (2.14c). Its most 



interesting feature is its frequency structure that fol- 

f2l (2) 

lows from {x ab '(x 1 , x 2 lx y cd '{x 3 , x 4 )} c dis being a function 



of t\ — t 2 and t 3 — T4O This implies that there are ef- 
fectively two free r integrals in $4 ^ , unlike the case of 



$4 , where there is only one. In frequency space, $4 



a (2) 



can be written 

1(2) 



>[M] 



niri2 



1 

8 

,(2) 



c?xi dx 2 dx 3 dx4 

(2), 



{Xab ( x i . x 2 ; w m ) X cd (x 3 , x 4 ; w„ 2 ) } d 



M° ( Xl , w ni ) (x 2 , -w ni ) Af c Q (xa , c^„ 2 ) 
xAf^(x 4 ,-w„ 2 ) . 

V«(l) 



(2.16) 



We see that $ 4 1 ; carries an extra factor of T compared 



to $ 4 ^ . The correlation function in Eq. ( 2 . 1 6j ) has,been 
calculated for noninteracting disordered electrons .O The 



important result of these authors was that the Fourier 
transform of the correlation function at u> ni = a>„ 2 = 
is finite. We have convinced ourselves that including 
interactions does no' 
then be replaced by 



interactions does not change this result. $4 ^ [M] can 



$:^[M] = -^|dx^|M«(x,c ril )| : 



nin 2 
X 



|M"(x, 



(2.17) 



with U4 a finite coefficient. The physical meaning of $4 
is easily determined. Consider the local in space and time 
contribution to the disorder average of the coefficient 
in Eq. ( ^8fj ), i.e. the term t = 1/T t ~ Xo(q = 0), 
which determines the distance from the critical point in 
the Gaussian theory. If we make to a random function of 
space with a Gaussian distribution and integrate out that 
randomness, then we obtain a term with the structure of 
<I>4 This term in our action thus represents a 'random 
mass' term, reflecting the fluctuations in the location of 
the critical point mentioned in the introduction. 

We conclude this section with a discussion of why 
the diagram shown in Fig. |l| gives the most important 
contribution to the correlation function {x }dis- The 
field theory of Ref. || allows for a systematic loop ex- 
pansion for these correlation functions, with the diagram 
shown in Fig. |l| the one-loop contribution. If we consider 
higher- loop corrections to this, we need to distinguish be- 
tween skeleton diagrams and insertions. Insertions will 
produce finite (in d > 2) renormalizations of the one- 
loop result. In the skeletons, each additional loop adds 
one frequency-momentum integral, and one independent 
propagator. Since,the propagators are at most a diffu- 
sion pole squared,B this means that each loop adds effec- 
tively a factor of q d ~ 2 to {x }dis> and hence we need 
to consider the lowest loop contribution for a given type 
of term. For the static terms, i.e. those of zeroth order 
in the frequency fl n , the lowest nonvanishing contribu- 
tion is at one-loop order since the zero-loop terms are 
of 0(f2 n ). For instance, the zero-loop contribution to 
{x^}dis is ~ |f2„|/(|f2„| +Dq 2 ), while the o ne-loop con- 
tribution adds the static piece to yield Eq. ( 2.11a ). We 
thus conclude that the most infrared divergent contribu- 
tion to the static Z-point correlation function is given by 
the diagram shown in Fig. [I] with dressed internal prop- 
agators, and for just getting the power of the leading 
divergence it suffices to use bare internal propagators. 
In the critical limit the zero-loop contr ibutio n to 



{x^}dis is proportional to Q n /q 2 , see Eq. ( 2.11 b| ) . In 



the next section we will discuss a fixed point where the 
dynamical exponent is z = d, i.e. fi n scales like q d at crit- 
ically. In the critical limit, the term of 0(fl n ) therefore 
show s the same scaling behavior as the term of 0(fi°), 
Eq. ( [2.12a[ ). The same is true for all higher correlation 
functions. It is easy to show that the zero-loop contribu- 
tion to {x^}dis goes like fl n /q 2l ~ 2 , while the one-loop 
contribution goes like q d ~( 2l ~ 2 K Anticipating again that 



G 



fi„ scales like q d at criticality, these two contributions 
have the same scaling behavior for all values of I. All 
higher loop contributions will be less important, as we 
have seen above. Also, higher orders in an expansion in 
powers of f2 can be neglected, since il n ~ q 2 in the dif- 
fusive propagators. Hence, each additional power of fi ra 
will lead to a factor that scales like q d ~ 2 at criticality 
and is less relevant than the terms of 0(f2°) and 0(O n ). 
We thus conclude that the most relevant term of 0(M l ) 
in the GLW functional has a coefficient ui that behaves 
effectively like 



ui 



d-2(l-l) 



(2.18) 



This covers both the terms of 0(^°) and 0(f2 n ), and 
all higher powers of £l n are less i mpor tant. The leading 
term for q — > 0, shown in Eq. ( 2.18| ), comes from the 
contribution of the l-th moment to the cumulant 
i.e. from {x }dis> while the other contributions, i.e. the 
subtraction terms in a given cumulant, are less divergent. 
However, due to the frequency structure of these terms 
that was explained above using the example of X^ 4 \ it 
is not obvious that they are unimportant for the critical 
behavior. We will analyze this point in the next section, 
and will find that only the subtraction term in t he q uar- 
tic cumulant, which is given explicitly in Eq. fl2.17|) , is 
important in that respect and thus needs to be kept. 



III. THE CRITICAL BEHAVIOR 

In the first part of this section we discuss the critical 
behavior of t he G aussian part of the theory defined by 
Eqs. (|2.8[ ) - (2.17). The renormalization group proper- 
ties of the Gaussian fixed point are also discussed. We 
then analyze the non-Gaussian terms in the field theory, 
and show that they are irrelevant, in the renormalization 
group sense, with respect to the Gaussian fixed point for 
all dimensions d > 2, except for a marginal operator in 
d = 4. This implies that the Gaussian theory yields the 
exact critical behavior for all of these dimensions, except 
for logarithmic corrections to scaling in d — 4 and d = 6 
that are discussed in Appendix B. We then construct the 
equation of state near the critical point, which requires 
a more detailed knowledge of the non-Gaussian terms in 
the field theory since it is determined in part by danger- 
ous irrelevant variables. We conclude this section with a 
discussion of the specific heat near the quantum critical 
point, and with a discussion of the quantum-to-classical 
crossover behavior. 



A. The Gaussian fixed point 



According to Eqs. (2.8) - (2.12), the Gaussian part of 
$ a is, 



d>«[M]=i££M"(q, w „) [t 



+ a d - 2 q 



d-2 



q u n 



+ a 2 q 2 + a w \uJ n \/q 2 } ■ M a (-q, -u n ) , (3.1s 



where 



to 



1-T< Xs (q^0,cj„ = 0) 



(3.1b) 



is the bare distance from the critical point, and ad-2, a 2l 
and duj are positive constants. 

W e first analyze the critical behavior implied by Eqs. 
( |3.l|) . Later we will show that for d > 2 fluctuations 
are irrelevant, and the critical behavior found this way is 
exact for these dimensions. By inspection of the Gaussian 
LGW functional in Eq. (3.1a) one obtains, 



l/(d-2) for2<d<4 
1/2 ford>4 



i-d 


-it 



for 2 < d < 4 
for d > 4 

for 2 < d < 4 
for d > 4 



(3.2a) 



(3.2b) 



(3.2c) 



Here v is the correlation length exponent, defined by 
£ ~ t~ v , with t the dimensionless distance from the 
critical point, r/ is the exponent that determines the 
wavenumber dependence of the order parameter suscep- 
tibility at criticality, (M Q (q, 0) M o (-q,0)> ~ q~ 2+v . z is 
the dynamical scaling exponent that characterizes critical 
slowing down by relating the divergence of the relaxation 
time, T r , to that of the correlation length, r r ~ £ z . 

Let us discuss, for later r efer ence, the critical behav- 
ior given by Eqs. (3^) and (3.2) from a renormalization 
group point of view. Let b be the renormalization group 
length rescaling factor. Under renormalization, all quan- 
tities change according to A — * A(b) — A, with [A] 
the scale dimension of A. The scale dimension of the 
order parameter is, 



[M(q,w n )] 



-1 + 77/2 



or, equivalcntly, 



[M(x,r)] = (d + 2)/2 



(3.3a) 



(3.3b) 



At the critical fixed point, a w and either ad-2 (for 
2 < d < 4), or ci2 (for d > 4) are not renormalized, 
i.e. there scale dimensions are zero. Using this, and 
[q] = 1, \c u n ] = z immediately yields Eqs. ( |3.2b| , 



3.2c) 



Equation (3.2a) follows from the relevance of £o> or its 
renormalized counterpart, t, at the critical fixed point. 
That is, the scale dimension of t is positive and given by 
\jv= [t] =2-r) = d-2. 



7 



B. The non-Gaussian terms 

We now show that all of the non-Gaussian terms in 
the field theory are renormalization group irrelevant with 
respect to the Gaussian fixed point discussed in the last 
subsection . Let us first consider the term of order M 3 . 
From Eq. ( fc.l3c|) we see that the most relevant coefficient 
at that order is ul d ^ for d < 4, and ul " 1 for d > 4. For 



simplicity, we will use the symbol U3 for the most relevant 
coefficient in a given dimension, i.e. U3 denotes u 



(d-4) 



for d < 4, and m 3 for d > 4. For d < 4 we then have 
to assign a scale dimension to the cutoff wavenumber q. 
The most obvious choice is to identify q with the inverse 
correlation length, which makes [q] = 1. While we 
will see later that this identification is not correct for all 
values of d, it will turn out that it provides an upper limit 
for [q]. Assu ming, then, [q] < 1, which will be justified 
in Sec. |EII C below, and using Eq. (3J3), we see that the 
scale dimension of the effective coefficient of the term of 
order M 3 is bounded by, 



[u a ] < -(d-2)/2 



(3.4) 



Similarly, using Eqs. ( 2.1 5d ) and (B.3), one sees that 
the coefficient 114 in Eq. ( p. 15a ) has a scale dimension 
bounded by, 



f-(d-2) for2<d<4 
M < < -2(d - 3) for 4 < d < 6 
I -d for d > 6 



(3.5a) 



and that the scale dimension of V4 in Eq. (2.17) is, 

[ V4 ] = -\4-d\ . (3.5b) 

Here a gain we denote the most relevant co effic ient in Eq . 
fl2.15d| ) by 1*4 for simplicity. Equations {3A) and ( |3.5| ) 
imply that the cubic and quartic terms in the field theory 
are renormalization group irrelevant with respect to the 
Gaussian fixed point for all d > 2 except for d = 4, 
where the coefficien t V4 is a marginal operator. As we 
explained after Eq. ( 2.17 ), V4 reflects the 'random mass' 
contribution to the action of the disordered magnet. We 
stress that this operator, while being strongly relevant 
with respect to the clean Gaussian fixed point, is not 
relevant with respect to the Gaussian fixed point of the 
present LGW functional. 

It is an easy matter to analyze the behavior of the 
higher order terms in the Landau expansion. The co- 
efficient of the general term of order M l in the LGW 
functional has a scale dimension bounded by, 



r -(d- 2) (I - 2)/2 for2<d<4 
M < \ -{l-2)-l(d-4)/2 for4<rf<6 
{-(l-4)-d(l-2)/2 ford>6 



(3.5c) 



This holds for the pure moment contribution to the Z-th. 
cumulant, and we see that all of the ui are irrelevant for 
d > 2. For the coefficients of the subtraction terms in the 



cumulants one easily convinces oneself that, while their 
scale dimension initially increases with d increasing from 
d = 2 (as does the scale dimension of V4, see Eq. (3.5b)), 
it stays negative for all d > 2. All of these terms are 
therefore irrelevant operators as well. We conclude that 
for d > 2 there is a critical Gaussian fixed point corre- 
sponding to a phase transition with Gaussian exponents 
(except for logarithmic corrections to the Gaussian crit- 
ical behavior in d = 4 due to the marginal operator v, 
and, as we will see later, also in d — 6), while for d < 2 
the non-Ga ussia n terms are potentially relevant. W e wil l 
see in Sec. [II C| below that the inequality in Eq. (3.5c) 
is actually an equality for 2 < d < 4, so that d — 2 is the 
upper critical dimension, and in d = 2 an infinite number 
of operators become marginal. 



C. The magnetization, and the magnetic 
susceptibility 

Although U4 is an iptelevant operator for d > 2, it is 
dangerously irrelevantliij for the magnetization, since it 
determines the critical behavior of the average magne- 
tization, to, as a function of t and an external magnetic 
field H . The technical reason for this is that to is a singu- 
lar function of u for u — ► 0. Schematically, the equation 
of state in mean-field theory is of the form, 



t TO + U4 TO = H 



(3.6) 



where we have s uppre ssed all numerical prefactors. Ac- 
cording to Eq. ( 2.15d ), 114 diverges for d < 6 as q d ~ 6 . 
Furthermore, if we were to keep higher terms of order m l 
in the equation of state, their coefficients would diverge 



as q 



d-2l 



This impli es t hat the cutoff q scales like m 



1/2 



and effectively Eq. (|3.6|) reads, 



t TO 



,(°), 



= H 



(3.7) 



with U4 from Eq. (^.15d ), and u another finite coef- 
ficient. From Eq. (3.7) we immediately obtain the ex- 
ponents (3 and 5, defined as m(t,H — 0) ~ t 13 , m(t = 

a, in - h 1 / 6 , 



f 2/(d-2) for 2 < d < 6 
\l/2 ford>6 



s= (d/2 for 2 < d < 6 
I 3 for d > 6 



(3.8a) 



(3.8b) 



In d = 6 logarithmic corrections to scaling occur, see 
Appendix B. 

Above we have used a general scaling argument to ob- 
tain Eq. (3/7). For small disorder, the same result can 
also be derived explicitly by mea ns of an infinite resum- 
mation. As mentioned after Eq. ( |3.6| ), the term of Q(m 3 ) 
with its divergent coefficient U4 ~ q d ~ 6 in Eq. ( |3.6| ) is 
only the first in an infinite series of terms that behave 
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like q' 



d-2l 



Calculating the prefactors of the divergent 



coefficients, one realizes that the divergencies are the con- 
sequence of an illegal expansion of an equation of state 
of the form 



t m 



const x mr 



V 



[(w„ + k 2 ) 2 + m 2 



H 



(3.6' 



Performing the integral one recovers Eq. (3.7). 

Next we determine the functional form of the equa- 
tion of state at nonzero temperature in order to ob- 
tain a scaling equation for to as a function of t, H, 
and T . This can m ost ea sily be done by utilizing Eq. 
( |2.12b| ). From Eq. fl2.12a|) it follows that the im-term 
in the equation of state has a correction of the form 
to (q 2 + T)(rf-2)/2 ^ m ( m + T)(d-2)/2 where we have 

used q 2 ~ to as ex plain ed after Eq. (3.6). Similarly, the 
term u m 3 in Eq. (|3.6|) with u ~ q d gets re plac ed by 



mac* 

(to + T)( d ~ 6 )/ 2 . At T = we recover Eq. (|J). For 
T < to in suitable units, there are corrections of 0(T/m) 
to the term ~ m d l 2 in that equation, while for T > to, i 
gets replaced by t + T l / 2u . All of these observations can 
be summarized in the following homogeneity law, 

m{t, T, H) = b- plv m(tb xlv \Tb'>'' v , Hb Sf}/u ) , (3.9a) 



with 



2v 



(3.9b) 



Similarly, the magnetic susceptibility, Xm> satisfies a ho- 
mogeneity law, 

Xm {t, T, H) = Wt v xmitb 1 /" \Tb^ v , Hb s ^/ V ) , (3.10a) 
with 



1 = P{5-\) = 1 



(3.10b) 



The Eqs. (3.9) and (3.1C) warrant some discussion. 
The scaling of T in these equations follows directly from 
Eq. ( |2.12b| ). The effective scale dimension of T in to and 
Xm is therefore 2 and not z. The salient point is that 
z is determined by the scaling of Q n or T with q in the 
Gaussian action, and hence in the critical propagator. 
However, the magnetization is calculated at tt n = q = 0, 
and its leading temperature dependence is determined by 
the diffusive modes, which feature fi„ ~ T ~ q 2 , rather 
than by the criti cal o nes. This leads to [T] — 4>/v, with <f> 
as given in Eq. ( 3.91: ). The proper interpretation of <f) is 
that of a crossover exponent associated with the crossover 
from the quantum to the thermal fixed point that occurs 
at any T > 0. Since z > <f)/u, the critical scaling T ~ b z 
would be the dominant temperature dependence if to and 
Xm depended on the critical modes. That they do not 
can also been seen from a determination of the magnetic 
susceptibility directly from the Gaussian ac tion: Recog- 
nizing that the coefficient of M 2 in $2, Eq. (3.1a), is the 
inverse spin susceptibility, and using Eq. (2.12t) again, 
we obtain 



Xm(t,T) = 



1 



t + TV 2 » 



(3.10c) 



in agreement with Eqs. (3.10a, 3.10b] ). 

We are now in a position to determine the exact scale 
dimension of U4, and of the other coefficients in the field 
theory, and to thus verify the assumption ma de i n the 
last subsection. As we have seen after Eq. ([3~6]), the 
cutoff q scales like to 1 / 2 , so that [q] = [m] e ff/2, where 
[m]eff = P/v is the effective scale dimension of to, i.e. 
the scale dimension after the effects of the dangerous ir- 
relevant variable U4 have been taken into account. From 



Eqs. ( pia[ ), ( pa] ) we see tha t \m] e g < 2, which justifies 



the assumption made in Sec. Ill B that led to the upper 
bounds on the scale dimensions of the ui . Repeating the 
power co unting argu ments that led to the inequalities, 
Eqs. (3.4) and ( 3.5a ), we obtain 



[u 3 ] = -(d-2)/2 , 



and 



f-(d-2) 
M = { ~{d 2 -l2)/(d- 
{ -d 



2) 



for 2 < d < 4 
for 4 < d < 6 
for d > 6 



(3.4') 



(3.5a') 



The exact scale dimension of the general coefficient ui , if 
desired, can be obtained by t he sa me argument. Again, 
the upper bound given in Eq. ( 3.5c ) is equal to the exact 
value for 2 < d < 4. These results complete the proof 
that d = 2 is the upper critical dimension for our prob- 
lem, and that an infinite number of operators become 
marginal in d = 2, and relevant in d < 2. 

Finally, we me ntion how to reconcile the value for the 
exponent (3, Eq. (3.8a), with scaling. Putting T = H = 
for simplicity, but keepi ng the dependence of to on U4 
explicitly, and using Eq. ( |3.3b ), we write 



m(t,u 4 ) = 6- (d+2)/2 TO(tfe 1 ^,u 4 6 [t141 ) 



Fro m Eq. (3.6) we know that m 



-1/2 



(3.9a') 
Using this in 

Eq. ( 3.9a ) changes the scale dimension of to from [to] = 
(d + 2)/2 to an effective value, [m] e // = [to] — [114] /2. 
With Eq. (3.5a') we then obtain [m] e // = P/u with the 
correct values of P and v. 



D. The specific heat 

The scaling equation for the specific heat is determined 
by the sum of the mean-field and the Gaussian fluctua- 
tion contribution to the free energy density, /. The mean- 
field contribution follows immediately from Eq. (|3.7|). 
The Gaussian fluctuation contribution, f@, which gives 
the leading nonanalyticity at the_critical point, can be 
calculated by standard methods.!!!] Neglecting an unin- 
teresting constant contribution to /q, we obtain, 
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2 In {H/m + aa-2 q 



d-2 



+a 2 q 2 + a u \u n \/q 2 ) 
+ In (x d H/m -{x d -l)t+ a d _ 2 q d ~ 2 

+a 2 q 2 + a u \uj n \/q 2 ) 



(3.11) 



Here x d = d/2 for 2 < d < 6, and x d = 3 for d > 6. The 
specific heat coefficient 7y is conventionally defined by, 



7v- = c v /T = -8 2 f/dT 2 



(3.12) 



Again we are interested only in scaling prop erties and 
not in exact coefficients. Schematically, Eqs. (3.11) and 
( HI give, 



■ /o 9 T+qd + q* + Hqym ' 



(3.13) 



with A an ultraviolet cutoff. 

Several points should be noted. Fi rst, f or all d that 
obey 2 < d < 4, j v g iven by Eq. (3.12), or exactly 



by Eqs. ( 3.11 ) and ( 3.12 ), is logarithmically singular for 
T, H — > 0. This d-independent logarithmic singularity is 
somewhat unusual. Wegner has discussed how logarith- 
mic corrections to scaling arise if arS«t of scale dimensions 
fulfills some resonance condition.E2l In the present case 
the appearance of a logarithm can be traced to the fact 
that the scale dimension of the free energy, d + z = 2z 
for 2 < d < 4, is a multiple of the scale dimension of 
T, which is z = d in this region. The unusual feature of 
the logarithm appearing in a range of dimensions, rather 
than only for a special value of d, is due to the dynam- 
ical exponent being exactly d in that range. Second, 
in Eq. (3.13) two different temperature scales appear. 
The first two terms in the denominator indicate that 
T ~ as one would expect from dynamical scaling. 
However, the last term in Eq. ( 3.13j ) contains the mag- 
netization, which in turn depend s on t he crossover tem- 
perature scale, T ~ £~ 2 , see Sec. HI C above. These two 



features imply that the scaling equation for 7^ should be 
written, 

<y v (t,T,H) = 9(4 -d) In 6 

+F 1 {tb 1/u ,Tb z ,Hb S(3/u ,Tb 2 ) . (3.14) 

Since z > 2, one can formally ignore the fourth entry in 
the scaling function since it is subleading compared to 
the second entry. The corrections to the resulting theory 
can be considered as 'corrections to scaling'. Notice that 
in contrast to the magnetization and the magnetic sus- 
ceptibility, the specific heat does depend on the critical 
modes, and hence contains the critical temperature scale. 
As mentioned in the last subsection, the latter is domi- 
nant when it is present, and jv provides an example for 
that. 



IV. DISCUSSION 

A. Theoretical aspects 

In this paper we have shown that in disordered itin- 
erant quantum ferromagnets the diffusive nature of the 
electrons leads to long-range interactions between spin 
fluctuations in an order parameter field theoty. As for 
classical models with long-range interactions ,EJ the crit- 
ical behavior of this field theory can be determined ex- 
actly. For d > 6, standard mean- field results are ob- 
tained, but for 2 < d < 6 one finds nontrivial critical 
behavior with dimensionality dependent exponents. For 
d < 2 our approach breaks down because the electrons 
are localized, and because in our field theory an infinite 
number of operators become marginal in d = 2. Th e ex - 
act c ritical expon ents for d > 2 are given by Eqs. fl3.2j), 
fl3.8|), and (3.10b), and the scaling properties of some of 
the more interesting physical properties are given by Eqs. 
d3.9[), (3.10a), and (3.14). In this subsection we discuss 
various aspects of these results that have not been cov- 
ered yet. 

First of all, there is an important conceptual question 
that should be considered. In our approach we have as- 
sumed that it is sensible to construct an order parameter 
field theory to describe the critical behavior of the order 
parameter. In general this procedure will break down 
if there are other soft or slow modes that couple to the 
order parameter fluctuations. That is, a more complete 
low energy theory might be needed. In the present case, 
the diffusons that lead to the n onan alytici tics in the bare 
field theory and to, e.g., Eqs. ( 2.12| ) and ( 2.15 ) are such 
soft modes. One should ask why we were able to pro- 
ceed with an order parameter description anyway, with- 
out running into unsurmountable difficulties due to the 
additional soft modes. A technical answer is that the dif- 
fusons did create problems, namely divergent coefficients 
in the LGW functional, but that for the present prob- 
lem these difficulties could be dealt with. Nevertheless, 
one might wonder what the theory would look like if the 
critical modes and soft diffusion modes were treated on 
a more equal footing. 

Without realizing it, we have previously addressed the 
above question. Early work on the metal-insulator tran- 
sition in disordered interacting electron systems showed 
that in two-dimensional systems without impurity spin- 
flip scattering, the triplet interaction scaled to large val- 
ues under renormalization group iterations.cj This was 
interpreted, incorrectly as it turned out later, as a_signa- 
ture of local moment formation in all dimensions .0 Sub- 
sequently, the present authors studied this problem in 
some detailo We were able to explicitly resum the per- 
turbation theory and show that at a critical value of the 
interaction strength, or of the disorder, there is a bulk, 
thermodynamic phase transition. The physical meaning 
of this phase transition was obscure at the time since no 
order parameter had been identified, and its description 
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was entirely in terms of soft diffusion modes. However, 
the critical exponents obtained are identical to those ob- 
tained here for the quantum ferromagnetic phase transi- 
tion, and in both cases logarithmic corrections to scaling 
are found. Because the exponents in the two cases are 
identical, we conclude that the transition found earlier 
by us, whose physical nature was unclear, is actually the 
ferromagnetic transition. One also concludes that our 
speculations about the nature of the ordered phase as 
an 'incompletely frozen spin phase' with no long-range 
magnetic order, were not correct. On the other hand, 
the techniques used in Ref. U allowed for a determina- 
tion of the qualitative phase diagram as a function of 
dimensionality, which our present analysis is not capable 
of. The theory given here not only explains the nature 
of the transition, but also explains why the critical be- 
havior at that phase transition could be obtained exactly 
in three dimensions: The long-range nature of the inter- 
actions between the order parameter fluctuations makes 
the critical phenomena problem exactly soluble. It is 
also interesting to note that the list of scaling scenar- 
ios for soft-mode field theories for disordered interacting 
fermions given in Sec. IV of Ref. || included the present 
case, namely a transition to a ferromagnetic state with 
an order parameter exponent (3 = 2v for d < 4. 

It should also be pointed out that our earlier theory 
depended crucially on there being electronic spin con- 
servation. This feature would be lost of there were some 
type of impurity spin flip scattering process. In that case, 
the soft modes that lead to the long-range order param- 
eter interactions acquire a mass or energy gap, and at 
sufficiently large scales the interactions are effectively of 
short range. The asymptotic critical phenomena in this 
case are described by a short-range, local order param- 
eter field theory with a random mass, or temperature, 
term. In this case the random mass term is relevant with 
respect to the Gaussian fixed point analogous to the one 
discussed here, which underscores the important role that 
is played by the order parameter being conserved in our 
model. The quantum phase transition in a model where 
it is not, and where the random mass term is therefore rel- 
evant with respect to the analogue of our Gaussian fixed 
point, is discussed elsewhere.t3 We also mention that the 
effect of fermionic soft modes on the ferromagnetic phase 
transition in clean systems has been discussed recently in 
Ref. §. 

Even though we quote exponent values formally for all 
d > 2, it should be mentioned that the region of valid- 
ity of our analysi s shr inks to zero as d approaches 2. As 
mentioned in Sec. II C , the reference system has all of the 
characteristics of the system described by the full action 
S 1 , except that it must not undergo a phase transition lest 
our separation of modes that is implicit in our singling 
out S* nt for the decoupling procedure break down. This 
requirement puts restrictions on the parameter values for 
which our procedure works. For instance, we cannot go 
to arbitrarily large disorder (at fixed d) without trigger- 
ing a metal-insulator transition within the reference sys- 



tem. For d — > 2 the metal-insulator transitions occurs at 
smaller and smaller values of the disorder, and in d = 2 
one obtains a very complicated, and unsolved, situation 
where various fluctuations compete with each other. 

We finally discuss why some of our results are in dis- 
agreement with Sachdev'scl recent general scaling anal- 
ysis of quantum phase transitions with conserved or- 
der parameters. For instance, it follows from our Eqs. 
( |3.1f4 |3T4| ) that the Wilson ratio, defined as W = 
(mjH)JJCv/T), diverges at criticality rather than be- 
ing a universal number as predicte d in R ef. |^. Also, for 
2 < d < 4 the function F 1 in Eq. ( |3.14| ), for t = and 
neglecting corrections to scaling, is a function of T/H, 
in agreement with Ref. ||, but for d > 4 this is not the 
case. The general reason for this breakdown of general 
scaling is that we work above an upper critical dimen- 
sionality, and hence dangerous irrelevant variables have 
to be considered very carefully, and on a case by case 
basis. This caveat is particularly relevant for quantum 
phase transitions since they tend to have a low upper 
critical dimension. It is well known that a given irrel- 
evant variable can be dangerous with respect to some 
observables but not with respect to others. Specifically, 
in our case the dangerously irrelevant variable U4 affects 
the leading scaling behavior of the magnetization, but 
not that of the specific heat coefficient, which leads to 
the divergence of the Wilson ratio. A simple example 
of this phenomenon is provided by classical <fi -theory in 
d > 4, where the dangerous irrelevant variable u (the 
coefficient of the </> 4 -term) affects the scale dimension of 
the magnetization, but not that of the specific heat.Ej In 
classical </> 4 -theory this point is obscured by the fact that 
the saddle-point contribution to the specific heat contains 
a discontinuity. This is often expressed as a — 0, with 
a the specific heat exponent. However, the approach to 
the discontinuity is described by the a suggested by hy- 
perscaling, namely a = 2 — d/2, see ch. VII. 4 of Ref. 
p^ . At the quantum FP the situation is clearer, since 
the saddle-point contribution to CV is subleading. It is 
also important to remember that different arguments of 
a scaling function can be affected in different ways by 
one and the same dangerous irrelevant variable. Here, 
the effective scale dimension of H in the specific heat is 
changed by 114 (from (3c? — 2)/2 to d in 2 < d < 4), but 
that of T is not, since U4 imports only the subleading 
crossover te mpera ture scale into jv via the appearance 
of m in Eq. (3.13). 



B. Experimental aspects 

In order to apply our theoretical results to experi- 
ments, one needs materials that show a transition from a 
paramagnetic state to a ferromagnetic one at zero tem- 
perature as a function of some parameter a;. Obvious can- 
didates are magnetic alloys of the form Fi_ x with P 
a paramagnetic metal, and F a ferromagnetic one. Such 
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FIG. 2. Schematic phase diagram for an alloy of the form 
PxFi- x . T c is the Curie temperature for the pure ferromag- 
net F, and x c is the critical concentration. 



materials show the desired transition as a function of the 
compositiori|-Barameter x; examples include F = Ni and 
P = Al, Ga.til The schematic phase diagram is shown in 
Fig. | 

One striking difference between our results for the 
quantum ferromagnetic phase transition at (x = x c , T — 
0) and the classical or thermal transition for Heisenberg 
ferromagnets are the numerical values of the exponents. 
For 3-d systems, our Eqs. (3.8) for instance predict /3 = 2, 
while the corresponding value for the thermal transition 
isf3 /3 c i ass ps 0.37. The large difference between the clas- 
sical and the quantum value should be easily observable 
if it was possible to measure the magnetization at a suf- 
ficiently low temperature as a function of x in order to 
observe the crossover between quantum and classical crit- 
ical behavior in the vicinity of x c . One possible way to 
do such an experiment would involve the preparation of 
many samples with different values of x over a small re- 
gion of x. It might also be possible to probe the magnetic 
phase transition by using the stress tuning technique that 
was used to study the metal-insulator transition in Si:P.E3 
Alternatively, one could prepare a sample with a value 
of x that is as close as possible to x c , and measure the 
magnetic field dependence of the magnetization, extrap- 
olated to T = 0, to obtain the exponent S. Again, there 
is a large difference between our prediction of 5 = 1.5 in 
d = 3, and the classical value <5 c i aS s ~ 4.86. 

Another possibility, that does not involve an extrap- 
olation to T = 0, is to measure the zero-field magnetic 
susceptibility as a function of both t = \x — x c \ and T. 
Equation (3.10a) predicts 



X, 



l (t,T)=T- 1 ' 2 f x (T/t 2 



(4.1) 



Here f x is a scaling function that has two branches, /+ 
for x > x c , and /~ for x < x c . Both branches ap- 




FIG. 3. Schematic prediction for a scaling plot of the mag- 
netic suscteptibility. 



proach a constant for large values of their argument, 
f x (y — > oo) = const. For small arguments, we have 
fx(y ~* 0) ~ Vv> while /~ diverges at a nonzero value 
y* of its argument that signalizes the classical transition, 

fx(v y*) ~ (y - y*)~ 7ci -, with 7class » 1.39 the 

susceptibility exponent for the classical transition. Our 
prediction is then that a plot of \ m T 1 / 2 versus T/t 2 will 
yield a universal function the shape of which is schemat- 
ically shown in Fig. |[ Notice that the exponents are 
known exactly, so the only adjustable parameter for plot- 
ting experimental data will be the position of the critical 
point. This is on sharp contrast to some other quantum 
phase transitions, especially metal- insulator transitions, 
where the exponent values are not even approximately 
known, which makes scaling plots almost meaningless. cj 
Finally, one can consider the low-temperature behavior 
of the specific heat. According to Eq. ( |3.14 ), as the tem- 
perature is lowered for x > x c the leading temperature 
dependence of the specific heat will be 



c v (T) ~TlnT 



(4.2a) 



At criticality this behavior will continue to T — 0, while 
for x > x r it will cross over to 



c y (T) - (kt) T 



(4.2b) 



For x < x c one will encounter the classical Heisenberg 
transition where the specific heat shows a finite cusp (i.e., 
the exponent a < 0). 
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APPENDIX A: THE WAVENUMBER 
DEPENDENT SPIN SUSCEPTIBILITY 

Here we calculate the wavenumber dependent spin sus- 
ceptibility in a disordered interacting Fermi system. In 
terms of the Q-matrix field theory reviewed in Ref. |[ the 
spin susceptibility reads, 



Xs(q, On) 



E 

r=0,3 



i-YT 



E 

mi ,7712 

r^mo — 



(q) 



X 



7722 — 71,7712 



(-q)> 



(Al) 



The classical matrix field Q comprises two fermionic de- 
grees of freedom, and a general matrix element JQ"f, has 
six indices: i is the spin index (with i — 1,2,3 the spin 
triplet channel), r = 0, 3 denotes the particle-hole chan- 
nel (r = 1, 2 would be the particle-particle channel), n, m 
are Matsubata frequency indices, and a, (3 are replica 
indices. The matrix elements with nm < describe the 
soft modes, i.e., the particle- hole excitations, while those 
with nm > are massive. For explicit calculations it is 
convenient to consider the nonlinear er-model version of 
the Q-field theoryj23'El where the massive modes are in- 
tegrated out. In the resulting effective model, Q 2 = 1, 
and a loop expansion can be set up by expanding all Q 
in terms of q n m = &(—nm)Q nm . The contribution of 
0(q 2 ) to Eq. ( |Al[ ) vanishes for fi„ = and hence does 
not contribute to Xo(q) = Xsi^-^n = 0). The one- loop 
contribution, which is the term of 0{q ), gives 



xo(q)~^I>I>»z>*(k)x> n (k) 



[2>„(k-q)-X£(k-q) 



(A2a) 



Here 2jVt and are the diffusive propagators of the 
theory.Ei Their structure is 



2? n (k) = l/(fc 2 + ^-^„) 



(A2b) 



with D a diffusion coefficient. has the same structure, 
with D replaced by D l ^ D. Since we are not interested 
in prefactors, we do not have to specify either D and 
or the prefactor i n Eq . (A2a). For the reasons discussed 
at the end of Sec. II C , the one-loop term suffices to cal- 



culate the leading infra red wavenumber dependence of 
Xo- Schematically, Eqs. jA^) yield, 



Xo(q) 



dk k 



d-l 



(A- 2 



(A3) 



with A an ultraviolet cutoff, from which one readily ob- 
tains Eq. fl2.12a| ). 



APPENDIX B: LOGARITHMIC CORRECTIONS 
TO SCALING IN D = 4 AND D = 6 

There are three distinct mechanisms that produce log- 
arithmic corrections to scaling: (1) Marginal operators, 
(2) Wegner resonance conditions between a set of scale 
dimensions, and (3) logarithmic corrections to the scale 
dimension of a dangerous irrclc>va,nt operator. The first 
two mechanism are well known.Ej The third is operative 
only above an upper critical dimension, and is therefore 
of particular interest for quantum phase transitions. 

In the present case, logarithmic corrections to scaling 
arise due to all three of these mechanisms. The second 
one produces corrections to the scaling of the specific 
heat in al l dimensions 2 < d < 4, as was discussed in 
Sec. [II D. The first one is operative in d = 4, where 
v 4 is a marginal operator, see Eq. ( p. 5b ). If desired, the 
resulting correctionS|-tp scaling can be worked out using 
standard techniques.! 2 !] Finally, the third mechanism pro- 
duces corrections to scaling in d — 6. According to Eq. 
(2.9), the coefficient 114 ~ Inq in d = 6. Via Eq. (3J3) or 
(3.7) this leads, for instance, to a leading behavior of the 
spontaneous magnetization, 



m(t,H = 0) 



t i/2 



and at the critical point we have 



m(t = 0,H) 



IH1/H)} 



1/3 



(Bl) 



(B2) 



Other consequences, e.g. for the specific heat in a mag- 
netic field, can be easily worked out. 
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